On the quasi-static evolution 
of nonequilibrium steady states 



Walid K. Abou Salem* 
Institute for Theoretical Physics, ETH Zurich 
CH-8093 Zurich, Switzerland^ 



Abstract 

The quasi-static evolution of steady states far from equi- 
librium is investigated from the point of view of quantum sta- 
tistical mechanics. As a concrete example of a thermodynamic 
system, a two-level quantum dot coupled to several reservoirs of 
free fermions at different temperatures is considered. A novel 
adiabatic theorem for unbounded and nonnormal generators of 
evolution is proven and applied to study the quasi-static evolu- 
tion of the nonequilibrium steady state (NESS) of the coupled 
system. 
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1 Introduction 

Recently, there has been substantial progress in understanding and 
rigorously proving the asymptotic convergence (as time t — ► oo) of a 
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state of a thermodynamic system, say one composed of a finitely ex- 
tended system coupled to one or more thermal reservoir, to a steady 
state, both in equilibrium ([JP1,2,BFS,FM,M1,2, DJ]) and far from 
equilibrium ([Rul,2,JP3,FMUe,MMSl,2]) from the point of view of 
quantum statistical mechanics. After the state of a certain thermody- 
namic system reaches a steady state, it is natural to ask how the state 
will evolve if the system is perturbed slowly over time scales that 
are large compared to a generic relaxation time of the system, and 
how much the state of the system will be close to the instantaneous 
(non)equilibrium steady state. 

This question was first addressed in [A-SF1], where the isother- 
mal theorem, an adiabatic theorem for states close to thermal equilib- 
rium, has been proven, and applications of this theorem to reversible 
isothermal processes have been discussed. Here, we pursue this ques- 
tion further by investigating the quasi-static evolution of states far 
from equilibrium from the point of view of quantum statistical me- 
chanics. 

According to the spectral approach to nonequilibrium steady states 
(NESS), the latter corresponds to a zero-energy resonance of the (ad- 
joint of the) C-Liouvillean; (see [JP3,MMS1,2]). Since the C-Liouvillean 
is generally nonnormal and unbounded, we prove an adiabatic theorem 
for generators of evolution that are not necessarily bounded or nor- 
mal. This theorem can be extended to study the adiabatic evolution 
of quantum resonances. [A-SF2] 

As a concrete example of a thermodynamic system, we consider 
a system composed of a two-level quantum system coupled to several 
fermionic reservoirs at different temperatures (for example, a quantum 
dot coupled to electrons in several metals). We apply the general adi- 
abatic theorem to study the adiabatic evolution of the NESS for this 
system. The main ingredients of our analysis are an adiabatic theo- 
rem for nonnormal and unbounded generators of evolution, a concrete 
representation of the fermionic reservoirs (Araki-Wyss representation 
[ArWy]), the spectral approach to NESS using C-Liouvilleans, and 
complex deformation techniques as developed in [HP,JP1,2,3]. 

The organization of this paper is as follows. In section 2, we state 
and prove a general adiabatic theorem (Theorem 2.2). This is the 
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key result of this section, which we apply in the subsequent sections 
to study the quasi-static evolution of nonequilibrium steady states. 
In section 3, we discuss the concrete physical model we consider: a 
two level quantum system coupled to several fermionic reservoirs at 
different temperatures. 1 In section 4, we study the C-Liouvillean corre- 
sponding to the coupled system using complex deformation techniques 
(Theorem 4.3), and recall the relationship between the NESS and a 
zero-energy resonance of the C-Liouvillean (Corollary 4.4). In section 
5, we apply Theorem 2.2 to study the adiabatic evolution of the NESS 
of the coupled system. The main result of this section is Theorem 
5.1. We also remark on the strict positivity of entropy production in 
the quasi-static evolution of NESS, and on a concrete example of the 
isothermal theorem [A-SF1]. Some technical details and proofs are 
collected in an Appendix. 
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2 A general adiabatic theorem 

So far, adiabatic theorems that are considered in the literature deal 
with generators of evolution which are self-adjoint; (see for example 
[AE]). This is expected, since the generator of dynamics in quantum 
mechanics, the Hamiltonian, is self-adjoint. However, for systems out 
of equilibrium, a generally nonnormal and unbounded operator, the so 
called C-Liouvillean, can be used to generate an equivalent dynamics 
on a suitable Banach space. Since we are interested in studying the 
quasi-static evolution of NESS, it is useful to prove an adiabatic theo- 
rem for nonnormal generators of time evolution. This is what is done 
in this section. 

lr Qie analysis can be directly generalized to the case when the small system is 
coupled to several bosonic reservoirs by using methods developed in [MMS1,2]. 
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Consider a family of closed operators {A(t)},t G R + , acting on a 
Hilbert space 7i. We make the following assumptions on A(t) in order 
to prove the existence of a time evolution and to prove an adiabatic 
theorem. All of these assumptions will be verified in the applications 
which are considered in the subsequent sections. 

(Al) A(t) is a generator of a contraction semi-group for all t G R + . 

(A2) A(t) have a common dense domain V C 7i for all t G R + . 

(A3) For z G p(A(£)), the resolvent set of A(t), let R(z,t) := (z - 
A(t)) -1 . Assume that R(—l,t) is bounded and differentiable as 
a bounded operator on 7i, and that A(t)R(— l,t) is bounded, 
where the Q stands for differentiation with respect to t. More- 
over, assume that for every e > 0, — e G p(A(t)). 

Let U (t) be the propagator that satisfies 

d t U(t)iP = -A(t)U(t)4> ,U(t = 0) = l, (1) 

for t > 0; ip G X>. We have the following result. 

Lemma 2.1 

Suppose that assumptions (A1)-(A3) hold. Then the propagator 
U{t) satisfying (1) exists and is unique, and \\U{t)ip\\ < ||-?/>||, for ip G 
V. 

The result of Lemma 2.1 is standard, and it follows from assump- 
tions (A1)-(A3) above and Theorem X.70 in [RS2]. 2 

2 Choose i] > and let U(t) be the propagator generated by A(t) := A(t) + r). It 
follows from (Al) that A(t) is a generator of a contraction semigroup. Furthermore, 
for t,t' € R + , A(t')A(t)~ 1 is bounded due to the closed graph theorem and (A2) 

(see [RSI]). Moreover, for small \t-t'\, ||(t'-t)(l(t')^W _1 -l)|| = \\A(t)A (t)\\ + 
o(\t - t'\), which is bounded due to (A3). By Theorem X.70 in [RS2] (or Theorem 
2, Chp XIV in [Yo], section 4), this implies, together with (Al) and (A2), that U(t) 
exists and is unique. In particular, ||J7(f)^|| < 1 uniformly in t > (for H^H = 1). 
We also have \\U(t)\\ = e^\\U{t)\\. Taking the limit r\ -» gives \U{t)i\)\ < 1. 
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Assume that A(t) = A(0) for t < 0, and that it is perturbed slowly 
over a time r such that A^ T \t) = A(s), where s := | G [0, 1] is the 
rescaled time. The following additional two assumptions are needed 
to prove an adiabatic theorem. 

(A4) The eigenvalue A(s) G (t(t4(s)) is isolated and simple, such that 

dist(\(s),a(A(s))\{\(s)}) > d, 

where d > is a constant independent of s G [0, 1], and A(s) is 
continuously differentiate in s G [0, 1]. 

(A5) The projection onto A(s), 



F\(s) ^- f U(z.s) ( l: . (2) 

'7a (s) 



2vri 



where 7a(s) is a contour enclosing A(s) only, is twice differen- 
tiable as a bounded operator. 

Note that, since A(s) is simple, the resolvent of A(s) in a neigh- 
borhood A/" of A(s) contained in a ball B(A(s),r) centered at A(s) with 
radius r < d is 

S) = + Ranalytic{z, s) , (3) 

where Ranaiytic{z, s) is analytic in A/". We recall some useful properties 
of the resolvent and the spectral projection P\(s); (see [Kal]). 

(i) It follows by direct application of the contour integration formula 
that 

(Px(s) f = P x (s) , (4) 

and hence 

P x (s)P x (s)P x (s) = . (5) 
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(ii) 

A(s)P x (s) = P x (s)A(s) = X(s)P x (s) . (6) 

Proof. 

A(s)P x (s) = -L / (A(s) - z + z)(z- A^y'dz 

= £ (s) + £ (s) ( 7^y + * 

= A(s)P A (s) , 

and similarly, P A (<s)^4(s) = A(s)P\(s)- 

(iii) It follows from (3) and (A4) that, for i] E C and | < |?7| < d, 
there exists a constant C < oo, independent of 77, such that 

||i2(A(a)+iM)|| <C, (7) 

uniformly in s G [0, 1]. Moreover, since (A(s) + rj) E p(A(s)), it 
follows by the spectral mapping theorem (see for example [Yo], 
Chp. VIII, section 7) and (A3) that R(\(s)+r], s) is differentiable 
as a bounded operator. 3 

We now discuss our general adiabatic theorem. Let U T (s,s') be 
the propagator satisfying 

d s U T (s, s') = -tA(s)U t (s, s') , U T (s, s) = l, (8) 

for s > s'. Moreover, define the generator of the adiabatic time evolu- 
tion, 

A a (s):=A(s)--[P x (s),P x (s)}, (9) 



3 We know that, for z,u € p(A), 

(z - A)- 1 = (1 + (z - w)(w - Ay 1 )- 1 ^ - A)- 1 . 

In particular, choose z = A(s) + 77 and u> = — 1. Differentiability of R(\(s) + rf) as 
a bounded operator follows from the latter identity and assumption (A3). 
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with the corresponding propagator U a (s, s') which satisfies 

d s U a (s, s') = -rA a (s)U a (s, s>) ; U a (s, s) = 1 , (10) 

for s > s'. 

By Lemma 2.1 and (A1)-(A3) and (A5), both propagators U T (s, s') 
and U a (s,s') exist and are unique, and ||C/ T (s, s')\\, \\U a (s,s')\\ < C for 
s > s', where C is a finite constant independent of s,s' G [0,1]. We 
are in a position to state our adiabatic theorem. 

Theorem 2.2 (A general adiabatic theorem) 

Assume (A1)-(A5). Then the following holds. 

(i) 

P x (s)U a (s,0) = U a (s,0)P x (0) , (11) 
for s > (the intertwining property). 

(ii) There is a finite constant C such that 

sup \\U T (s,0)-U a (s,0)\\<TJ-, 
se[o,i] 1 + r 

for t > 0. In particular, 

sup \\U T (s,0)-U a (s,0)\\=O(r- 1 ), 

se[o,i] 

as r — > oo. 

Proo/. 

(i) Equality holds trivially for s = 0, since C/ a (s, s) = 1. Let 

h(s,s') :=U a (s,s')P x (s')U a (s',0) , (12) 

for < s' < s. 
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Using (6), (10) , the definition of A a (s) and the fact that P x (s)P x (s)+ 
P x (s)P x (s) = P x (s), it follows that 

d s ,h(s,s') = d s ,(U a (s, s')P x (s')U a (s',0)) 

= TU a ( S ,s'){A a (s')Px(s')-Px(s')A a (s')}U a (s',0) 

+ U a (s,s')P x (s')U a (s',0) 
= U a (s,s'){-P x (s')P x (s>) 
-P x (s')P x (s') + P x (s')}U a (s',0) 
= . 

Therefore, 

h(s, s') = h(s). 

In particular, 

h(s, s) = h(s, 0), 

which implies claim (i). 

(ii) Consider ip e V, where the dense domain V appears in assump- 
tion (A2). We are interested in estimating the norm of the dif- 
ference (U T (s, 0) — U a (s, 0))^ as r — > oo. Using (8), (10) and the 
Duhamel formula, 

(U T (s,0)-U a (s,0))1> = - f ds'd s iU T (s,s')U a ( S ',0))4> (13) 

Jo 



JO 



ds'(U T (s, s')[P x (s'), P x (s')}U a (s', 0))^ 

(14) 

Let 

X(s) := -L / dzR(z, s)P x (s)R(z, s) , (15) 

where 7a(s) is a contour of radius centered at A(s), and 
where d appears in (A4). Then 

[X(s),A(s)} = ^-l dz[z-A(s),R(z,s)P x (s)R(z,s)} 

= P x (s)P x (s) - P x (s)P x (s) = [P x (s), P x (s)ll6) 
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Assumptions (A3),(A4) and the spectral mapping theorem im- 
ply that, for z G 7a(s) C p(A(s)), R(z, s) is different iable as a 
bounded operator. Together with (A5), this implies that, 

\\X(s)\\ < d, (17) 
\\X(s)\\ < C 2 , (18) 

where C\ and C 2 are finite constants independent of s G [0, 1]. 
Moreover, 

U T (s,s')[X(s'),A(s')]U a (s',0) = ±{-d s ,U T (s : s')X(s')U a (s',0) 

T 

+ U T (s, s')(X(s') [P x (s>), P x (s>)])U a (s>, 0) + U T (s, s')X(s')U a (s', 0)} 

Together with (16), one may write the integrand in (13) as a 
total derivative plus a remainder term. Using the fact that V is 
dense in H and (A5), 

||C/ T ( S ,0)-C/ a ( S ,0)||<- sup[^||X( S )||+^||X( S )||], (19) 

T »e[o,i] 

where C[,i = 1,2 are finite constants independent of s G [0, 1]. 
Together with (17) and (18), this implies 

sup ||[/ r ( s ,0)-C/ a ( S ,0)||<-^, (20) 
se[o,i] 1 + r 

for t > 0, where C is a finite positive constant. □ 



Next, we discuss a concrete model of a thermodynamic system to 
be studied subsequently. 

3 The Model 

As an example, we consider a two-level quantum system E coupled to n 
reservoirs , 7Z\, ■ ■ ■ , lZ n , n > 2, of free fermions in thermal equilibrium 
at inverse temperatures Pi, - ■ ■ , (3 n . 
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The small system 

The kinematical algebra of S is (9 s = Ai(C 2 ), the algebra of com- 
plex 2x2 matrices over the Hilbert space Tip = C 2 . Its Hamiltonian 
is given by if s = WQa^, where <Ji,i = 1,2,3, are the Pauli matrices. 
When the system E is not coupled to the reservoirs, its dynamics in 
the Heisenberg picture is given by 

4*(a) := e ^ E (*-*) ae -^ s (*-*) , (21) 

for a G s . 

A physical state of the small system, u; E , is described by a density 
matrix p s . We assume that ps > 0, ie, uP is faithful. The operator 
k s = p 1 ^ 2 belongs to the space of Hilbert-Schmidt operators, which 
is isomorphic to Hp <g> Hp . Two commuting representations of (9 s on 
lip <g> lip are given by 

7r s (a) := a® I s , (22) 
7r#( a ) : = I s ® C s aC s , (23) 

where C s is an antiunitary involution on Hp corresponding to complex 
conjugation in the basis of the eigenvectors of iJ s ; (see for example 
[BFS]). 

The generator of the free dynamics on the Hilbert space Hp ® VP 
is the standard Liouvillean 

Cp = <g> I s - I s <g> H s . (24) 

The spectrum of £ s is 0"(£ s ) = {— 2cuq, 0, 2o; }, with double de- 
generacy at zero. 

Let a> s be the initial state of the small system S, with correspond- 
ing vector f2 s e 7i s ® 7i s . The modular operator associated with up 

is A s = lu s ® cj s , and the modular conjugation operator, J s , is 
given by 

for E TiP. If cj s corresponds to the trace state, then A s = 1 S ®1 E . 
The reservoirs 
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Each thermal reservoir is formed of free fermions. It is infinitely 
extended and dispersive. We assume that the Hilbert space of a single 
fermion is h = L 2 (R + ; £>), where B is an auxiliary Hilbert space, and 
m{u)du is a measure on R + . We also assume that the single- fermion 
Hamiltonian, h, corresponds to the operator of multiplication by u G 
R + . For instance, for reservoirs formed of nonrelativistic fermions in 
R 3 , the auxiliary Hilbert space B is L 2 (S 2 ,da), where S 2 is the unit 
sphere in H 3 ,dcr is the uniform measure on S 2 , and u = \k\ 2 , where 
k G R 3 is the particle's momentum. In the latter case, the measure on 
R + is choosen to be m{u)du = \-Judu. 

Let b and b* be the annihilation-and creation operators on the 
Fermionic Fock space JF(L 2 (R + ; £>)). They satisfy the usual canonical 
anticommutation relation (CAR) 

{&#(/), b*(g)} = 0, (25) 
{b(f),b*(g)} = (f,g)l, (26) 

where b* stands for b or b*, /, g G L 2 (R + ;£>), and (•, •) denotes the 
scalar product in L 2 (R + ;£>). Moreover, let £l n denote the vacuum 
state in F(L 2 (R+; B)). 

The kinematical algebra, O nx , of the i th reservoir TZi, % — 1, • • • , n, 
is generated by bf and the identity l n \ The free dynamics of each 
reservoir (before the systems are coupled) is given by 

<(bf(f))^bf(e^- s >f), (27) 

fori = 1,--- G L 2 (R + ;B). 

The (a-jii, A)-KMS state, oj ni , of each reservoir TZi,i = 1, • • • ,n, 
at inverse temperature is the gauge invariant, quasi-free state uniquely 
determined by the two-point function 

<** i (W)b i (f)) = (f, Pfk {.)f), (28) 

where pp.(u) := 

Next, we introduce Ff w := ^(L 2 (R+; B)) ® ^(L 2 (R+; B)), 
the GNS Hilbert space for the Araki-Wyss representation of each 
fermionic reservoir TZi associated with the state uj n \ [ArWy]. Denote 
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by the vacuum state in J^ i (L 2 (R + ; B)), with = 0. The 

Araki-Wyss representation, 7Tj, of the kinematical algebra O ni ,i = 
1, • • • , n, on is given by 



where iVj = rfrj(l) is the particle number operator for reservoir IZi. 
Furthermore, Q Ui <g> Vt ni G JF/ 11 ^ corresponds to the equilibrium KMS 
state u ni of reservoir 7^.j. 

The free dynamics on the GNS Hilbert space J-f w of each reser- 
voir IZi is generated by the standard Liouvillean C n \ The modular 
operator associated with (O n \uj ni ) is given by 



and the modular conjugation is given by 

<g> $) = (-l)^- 1 )/ 2 ^ <g> (_i)^W-i)/2^ 

for $ G J^/ 11 ^; (see, for example, [BR]). 

In order to apply the complex translation method developed in 
[HP,JP1,2,3], we map T* w = ^(L 2 (R+; B)) ® ^(L 2 (R+; B)) to 
F ni (L 2 (Tl; B)) as done in [JP3]; (using the isomorphism between L 2 (R + ; £>)© 
L 2 (R + ; B) and L 2 (R; B), the latter having measure du on R). To ev- 
ery / G L 2 (R + ;£>), we associate functions fp,ff G L 2 (R;£>), with 
measure du on R, by setting 



h(VT^ f) ® i ni + (-1)* ® 6;(VPft /) 



(29) 



A 7 ^ 





(30) 



and 



\ y/m(-u)iy/l - pp(-u) f(-u, a),u<0 
= ifp(-u,<r), 
where m(u)du is the measure on R + , see eq. (29). 




(31) 



12 



Let a,i and a* be the annihilation and creation operators on T Ui (L 2 (R 
Then 

MUf) + W)) - Oi(/ A ) + (32) 
(W) + W)) - + (33) 

(g> n n \ (34) 

where fi 7 ^ is the vacuum state in ^ r7 ^ i (L 2 (R; B)). 4 

Moreover, the free Liouvillean on JF^(L 2 (R; 6)) for the reservoir 
IZi is mapped to 

L Ui = dTi{u) , (35) 

where «GR. 

The coupled system 

The kinematical algebra of the total system, S V Hi V • • • V 7£ n , 
is given by 

O = s ® O ni <g> • • • <g> O 7 ^ , (36) 

and the Heisenberg-picture dynamics of the uncoupled system is given 
by 

«o = a s <g> a^i ® • • • <S> . (37) 

The representation of O on W := <g> T~F <g> ^(L 2 (R; #)) <g> 
• • • <g> JF^"(L 2 (R; £>)), determined by the initial state 

u = u; s <g> cj^ 1 <g> • • • <g> (38) 

by the GNS construction, is given by 

7T = 7T S <g> 7Ti ® • • • <g) 7r n , (39) 

and an anti-representation commuting with 7r by 

7T # = 7rf <g) 7rf (g) • • • <g) 7T# . (40) 
4 For a discussion of this map, see Theorem 3.3 in [JP3]; (see also the Appendix). 
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Moreover, let Q := f2 s <8> Q Ul ® • • • <g) ft 11 " denote the vector in H 
corresponding to the state uj. Denote the double commutant of tt(0) by 
M. := tt(0)", which is the smallest von Neumann algebra containing 
n(0). 

For a G O, we abbreviate 7r(a) by a whenever there is no danger 
of confusion. The modular operator of the total system is 

A = A s <g> A^ 1 <g> • • • <g> A Un , 

and the modular conjugation is 

J = J s ® J w i <g> ...(g, J^». 

According to Tomita-Takesaki theory, 

J.MJ = M' , A tt A<A- ft = M , 
for £ G R; (see for example [BR]). Furthermore, for a G M , 

JA 1/2 afi = a*fi. (41) 
The Liouvillean of the total uncoupled system is given by 

n 

£ = £Z + J2£ n >. (42) 
i=i 

This defines a selfadjoint operator on Ti. 

The system E is coupled to the reservoirs 72. i, • • • , lZ n , through 
an interaction gV(t), where V(t) G O is given by 

n 

V(t) = J^Wi ® [W*(*)) + &:(/*(*))]} , (43) 

i=l 

(Ti, % — 1, 2, 3, are the Pauli matrices, and /; G L 2 (R + ; £>), « = 1, • • • , n, 
are the form factors. 

The standard Liouvillean of the interacting system acting on the 
GNS Hilbert space TC is given by 

C g (t)=Co + gI(t), (44) 
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where the unperturbed Liouvillean is defined in (42), and the interac- 
tion Liouvillean determined by the operator V(t) is given by 



I(t) = {V(t) - JV(t)J} 

n 

= ® i s ® KihAt)) + atifijkit))] 

- a s ® ax ® + a,(/i(0)]}, (45) 

where aj, a* are the annihilation and creation operators on the fermionic 
Fock space JF^(L 2 (R; B)). Note that since the perturbation is bounded, 
the domain of C g (t) is V(C g (t)) = V(C ). 

Let U g be the propagator generated by the standard Liouvillean. 
It satisfies 

d t U g (t,t') = -ijC g (t)U g (t,t') ;U g (t,t) = 1 , (46) 
for t > t' . The Heisenberg-picture evolution is given by 

af (a) = U* g (t, t')aU g (t,t') , (47) 

for a eO. 

Generally, the kernel of £ g (t),Ker £ g , is expected to be empty 
when at least two of the reservoirs have different temperatures. 5 . This 
motivates introducing the so called C-Liouvillean, L g , which generates 
an equivalent dynamics on a suitable Banach space contained in 7i 
(isomorphic to O) and which, by construction, has a non-trivial kernel. 

Consider the Banach space 

C(0,n) := {att:ae £>}, 

with norm || afi | |oo = ||a||. Since Q is separating for O, the norm \\afl 
is well-defined, and since f2 is cyclic for O, C(0, Q) is dense in H. 

5 This is consistent with the fact that the coupled system is not expected to 
possess the property of return to equilibrium if the reservoirs have different tem- 
peratures (or chemical potentials). One can verify that, indeed, this is the case 
when assumptions (Bl) and (B2), below, are satisfied; (see [JP3,MMS1,2]). 
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Let U g (t,t') be the propagator given by 

a t /(a)n = U g (t,lf)an, (48) 

and 

U g (t',t)Q = n. (49) 
Moreover, let L g (t) be its generator, ie, 

d t Ug(t, t') = iUg(t, t')Lg(t) With Ug(t, t) = 1. (50) 

Differentiating (48) with respect to t, setting t = t', and using 
(50), (47) and (41), one obtains 

[(Co + gV(t))a - a(C + gV(t))]Q = [(C + gV(t))a - (V(t)a*)*]Q 

= (C + gV(t) - gJA l l 2 V(t)/\- l l 2 J)att 
= L g (t)aVL . 

Hence, the C-Liouvillean is given by 

L g (t) := Co + 9 V(t) - gj A 1 / 2 V(t)A~ 1/2 J . (51) 

Note that, by construction, 

L g (t)Q = 0, 

for all t e R. 

Next, we discuss the assumptions on the interaction. For 5 > 0, 
we define the strips in the complex plane 

1(5) :={zeC: \Imz\ < 5} 

and 

I~(S) ■= {z E C : -5 < Imz < 0}. (52) 

Moreover, for every function / G L 2 (R + ; B) , we define a function / 
by setting 

16 



where m(u)du is the measure on R + . Denote by H 2 (5,B) the Hardy 
class of analytic functions 

h : 1(5) -> B, 

with 



H 2 (5,t3) '■= sup / \\h(u + i9)\\ 2 B du < oo. 

|0|<<5 JR 



We require the following basic assumptions on the interaction term. 

(Bl) Fermi Golden Rule. 
Assume that 

n 

^||^(2^,t)|| e >0, (54) 
i=\ 

for almost all t G R, which is another way of saying that the 
small system is coupled to at least one reservoir, to second order 
in perturbation theory. 

(B2) Regularity of the form factors. 

Assume that 35 > 0, independent of t and i — 1, • • • ,n, such 
that 

e -W 2 /,(«,t)e^S), (55) 

the Hardy class of analytic functions. This implies that the map- 
ping 

R 3 r -> A ir V(t)A- ir G M , (56) 

(where A = A s ® A^ 1 <g>- • -®A n " is the modular operator of the 
coupled system, and Ai = tt(0)",) has an analytic continuation 
to the strip 1(1/2) = {z G C : \Imz\ < 1/2}, which is bounded 
and continuous on its closure, Vt G R. 

(B3) Adiabatic evolution. 

The perturbation is constant for t < 0, V(t) = V(0), and then 
slowly changes over a time interval r such that V r (t) = V(s), 
where s = t/r G [0, 1] is the rescaled time. We also assume that 
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V(s) is twice differentiable in s G [0, 1] as a bounded operator, 
such that 

R9r^ A tr dlV(s)A- tr e M ,j = 0,1,2, (57) 

has an analytic continuation to the strip {z e C : |Jmz| < 1/2}, 
which is bounded and continuous on its closure. This follows if 
we assume that there exists 5 > 0, independent of s and % = 
1, • • • ,n, such that 

e-^ 2 dif l (u,s)eH 2 (5,B) , (58) 

the Hardy class of analytic functions, for j = 0,1,2. This as- 
sumption is needed to prove an adiabatic theorem for states 
close to NESS. 6 

Let U g be the propagator generated by the adjoint of the C- 
Liouvillean, ie, 

d t U g (t,t') = -iL* g (t)U g (t,t'), (59) 
U g (t,t) = 1. (60) 

Assumption (B2) implies that the perturbation is bounded, and hence 
the domain of L#, where L* stands for L g or L*, is 

V(Lf) = V(£ ), 

and U g , U g are bounded and strongly continuous in t and t' . 

4 The C-Liouvillean and NESS 

In [JP3, MMS1,2], it is shown that, when the perturbation is time- 
independent, and under reasonable regularity assumptions on the form 

6 When the reservoirs are formed of nonrelativistic fermions in R 3 , an example 
of a form factor satisfying assumptions (B1)-(B3) is given by 

f i (u,s) = h i (s)\u\ 1 ^e-W 2 , 
where is twice differentiable in s. 
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factors, the state of the coupled system converges to a nonequilibrium 
steady state (NESS) which is related to a zero-energy resonance of the 
adjoint of the C-Liouvillean. Here, we study the C-Liouvillean in the 
time- dependent case, and relate a zero-energy resonance to the instan- 
taneous NESS. The statements made in this section have been proven 
in [JP3] (see also [JP1,2]) for the time-independent case. Extending 
those results to the time- dependent case is straightforward, since we 
study the spectrum of the Liouvillean at each fixed moment of time. 
However, a sketch of the proofs of all the statements made in this sec- 
tion is given in the Appendix to make the presentation self-contained. 

We first study the spectrum of L* using complex spectral defor- 
mation techniques as developed in [HP,JP1,2,3]. 

Let Uj be the unitary transformation generating translations in 
energy for the i th reservoir, i = 1, • • • ,n. More precisely, for /j G 
L 2 (R," B) 

u i (0)/ i (u)=/f(u) = / i (u + 0). 

Moreover, let 

Ui{9) := T t (u, t (9)) 

denote the second quantization of Uj(0). 

Explicitly, Ui(9) = e~ tdAi , where Ai := idTi(d Ui ) is the second 
quantization of the generator of energy translations for the i t/l reservoir, 
% — 1, • • • , n. We set 

U(0) := I s <g> 1 E <g> U^e) (8) ••• (8) U n (6) . (61) 

Define 

L* g (t,6):=U(9)L* g (t)U(-9), (62) 

which is given by 

L;(t,9) = £ + N9 + gV tot (t,9) , (63) 

C = £ E + X\£^, C n > = dT(ui),i = 1, • • • , n, N = £\ N t , the total 
particle number operator, and 

v to % 9) = ® i E ® [a*(/fS(o) + - ^ ® (p s ) _1/2 

^(P S ) 1/2 (-l) JV 1a i (e ft ^ 2 /?(0) + <(e-^ /2 /5 e) (0)]} • 
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It follows from assuption (B2) that, for 9 G 1(5), V g tot (t, 9) is a bounded 
operator. Hence L*(t, 9) is well-defined and closed on the domain V := 
V(N)r)V(C Ul )r)- ■ ■rXD(C 1ln ). When the coupling g = 0, the pure point 
spectrum of £ is Cpp(£o) — {— 2cj , 0, 2w }, with double degeneracy 
at 0, and the continuous spectrum of C is cr cont (C ) = R. Let 

C (9) := Co + N9. 

We have the following two lemmas. 

Lemma 4.1 For 9 E C, the following holds. 

(i) For any ip G V, one has 

||£o(^l| 2 = ||£o(i?eW| 2 + \Im9\ 2 \\Ni/j\\ 2 . (64) 

(ii) If ImO 7^ ; then C (9) is a normal operator satisfying 

Co(9)* = C (9) , (65) 

andV(C (9)) = V. 
(Hi) The spectrum of C (9) is 

Vcont(£o(6)) = {n9 + s:ne N\{0} and s G R}, (66) 
a pp (C (9)) = {E r .j = 0,--- ,3}, (67) 

where E 01 — 0, E 2 — —2u and E 3 = 2u , (the eigenvalues of 
Lemma 4.2 

Suppose assumptions (Bl) and (B2) hold, and assume that (g, 9) G 
C x I~(5). Then, for each fixed time t G R, the following holds. 

(i) V(L* g (t,9)) = V and (L* g (t,9))* = L w (t,9). 

(ii) The map (g,9) — > L*(t,9) from C x I~(5) to the set of closed 
operators on H is an analytic family (of type A) in each variable 
separately; (see [Kal], chapter V, section 3.2). 
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(Hi) For j£R finite and Imz large enough, 

s - lim (L* g (t, 9) - z)- 1 = (L* g (t, Re9) - z)' 1 . (68) 



We now apply degenerate perturbation theory, as developed in 
[HP], to compute the spectrum of L*(t,9). Using contour integra- 
tion, one may define the projection onto the perturbed eigenstates 
of L* g (t,9), for 9 E I' (5). Let 

P 9 (t,0)^f^(z-L;(t,9))-\ (69) 

where 7 is a contour that encloses the eigenvalues Ej, j = 0, • • • , 3, at 
a distance d > 0, such that, for sufficiently small \g\ (to be specified 
below) the contour also encloses Ej(g,t), the isolated eigenvalues of 
L*(M).Welet 

p = i s ® i s ® \n n " (g)---(g)tt ni }(n ni <g> ...<g> 

where I s corresponds to the identity on 7Y S and VL ni corresponds to 
the vacuum state in JF^ 1 (L 2 (R; £>)). Furthermore, we define 

T g (t) :=P P g (t,9)P . (70) 

Consider the isomorphism 

S g (t, 9) := T-^ 2 (t)P P g (t, 9) : Ran(P g (t, 9)) -> Ran(P ) (71) 
and its inverse 7 

S;\t, 9) := P g (t, 9)P T- 1 / 2 (t) : Ran(P ) -> Ran(P g (t, 9)). (72) 
We set 

M g (t) := P P g (t, 9)L* g (t, 9)P g (t, 9)P , (73) 

7 It follows from (78), Theorem 4.3 (i) below, that T g (t) -> 1 on Ran(P g (t,6)) 
as 5 — > 0, and hence S fl (i, 0) is a well-defined operator on Ran(P g (t, 9)). By (70), 
it has the right inverse Sg 1 ^, 9). Moreover, dimRan(P g (t,6)) = dimRan(Po) for 
g small enough, and hence <5~ 1 (t, 0) is the inverse of S g (t, 9). 
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and define the quasi-C-Liouvillean by 

E p (f) := S g (t,9)P g (t,9)L;(t,9)P g (t,9)S; 1 (t,9) = T g ^(t)M g {t)T g ^(t). 

(74) 

Let k = min{5, ■ ■ • , ■£-}, where 5 appears in assumption (B2), sec- 
tion 3, and /3i, • • • , /3 n , are the inverse temperatures of the reservoirs 
IZi,--- ,1Z n , respectively. For 9 G I~(h) (see (52)), we choose a pa- 
rameter v such that 

k + I v\ 

-k < v < and - k < Im9 < ^— ( 75 ) 

We also choose a constant g± > such that 

<7iC < (fc - M)/2, (76) 

where 

C := sup \\V tot {t,9)\\ (77) 
6»e/(5),teR 

< sup ^ £ |1 + e-^l-V^sll^^H^^ + He-A^H^fcB)), 

t6R,z6/(<5) ^ { 

which is finite due to assumption (B2). 
Theorem 4.3 

Suppose that assumptions (Bl) and (B2) hold. Then, for g\ > 
satisfying (76), 9 G I~(k), v satisfying (75), and t G R /ixed, £/ie 
following holds uniformly in t, ie, g\ is independent of t. 

(i) If\g\ < gi, the essential spectrum of the operator L*(t, 9) is con- 
tained in the half-plane C\S(z/), where S(z/) := {z G C : Imz > 
u}. Moreover, the discrete spectrum ofL* (t,9) is independent of 
9 G I~(k). If \g\ < \g\, then the spectral projections P g (t,9), 
associated to the spectrum of L*(t,9) in the half-plane S(z/), are 
analytic in g and satisfy the estimate 

\\Pg(t,e)-P \\<l. (78) 
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(ii) If \g\ < y> then the quasi- C-Liouvillean £ g (£) defined in (74) 
depends analytically on g, and has a Taylor expansion 

CO 

Z g (t)=jf + J29 2j X {2j \t). (79) 
The first non-trivial coefficient in (79) is 

s(2)(t) = \ i t){z - £S) " + {z - c * r1 ^ *)) ' 

where £(z, t) := P V tot (t, 9){z - C (9)y 1 V tot (t, 9)P . 

In fact, one may apply second order perturbation theory to cal- 
culate the perturbed eigenvalues of L*(t,9). To second order in the 
coupling g, 

E (g,t)=0, 

Ek{g,t) = -z7rg 2 J2U( 2 ^t)\\l + 0(g 4 ) , 

i 

and 

E 2 ^t)= T {2^-\g 2 VV f ^^—^11^,0111) 
- l ^J2\m^,t)\\l + 0(g'), 

i 

where W denotes the Cauchy principal value (see the Appendix). 

The following corollary follows for the case of time-independent 
interactions; (see [JP3,MMS1,2]). 

Define 

D := I s <g> I s <g> e ^ kA ^i ® . . . (g, e -^„ ; ( 80 ) 

where := dT(y/p? + 1), and := i9 Ui is the generator of energy 
translations for IZi, i — 1, • • • ,n. Note that D is a positive bounded 
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operator on Ti such that Ran(D) is dense in 7i and DQ = Q. This 
operator will act as a regulator which is used to apply complex defor- 
mation techniques. Let a g = a*' . 

Corollary 4.4 (NESS) 

Suppose assumptions (Bl) and (B2) hold, and that the perturba- 
tion V{t) =V is time-independent. Then there exists g± > such that, 
for < \g\ < gi and afl G V(D^ 1 ), the following limit exists, 

lim (n, aUa)n) = (Q g , D^aQ) , (81) 

where Q g corresponds to the zero-energy resonance of L*, and a g is 
the perturbed dynamics. For a G O test , a dense subset of O (that will 
be specified below), this limit is exponentially fast, with relaxation time 
T R = 0(g-*) s 

Moreover, [JP3,MMS1,2] prove strict positivity of entropy pro- 
duction in the latter case, which is consistent with Clausius' formu- 
lation of the second law of thermodynamics. See [FMUe] for another 
proof using scattering theory of the convergence to a NESS and strict 
positivity of entropy production when two free fermionic reservoirs at 
different temperatures or chemical potentials are coupled. 

5 Quasi-static evolution of NESS 

In this section, we apply Theorem 2.2, section 2, to investigate the 
quasi-static evolution of NESS of the model system introduced in sec- 
tion 3. 

Together with assumption (Bl), we assume (B3), ie, V T (t) = 
V(s), where s G [0, 1] is the rescaled time with sufficient smoothness 
properties of the interaction. From Theorem 4.3, section 4, we know 
the spectrum of the deformed adjoint of the C-Liouvillean, L*(t, 9) = 
U(6)L* g (t)U(-6), for 9 G I'(k), where k = min(<5, and 5 

8 In fact, by assuming additional analyticity of the interacting Hamiltonian, one 
may show that this result holds for any initial state normal to u>; see [JP3,MMS1,2]. 
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appears in assumption (B3). Let 70 be a contour enclosing only the 
zero eigenvalue of L*(s, 9), for all s G [0, 1], and 

P ° M):= i ^-W'W 1 ' (82) 

the spectral projection onto the state corresponding to the zero eigen- 
value of L*(s, 9). Moreover, let h test = V(e k ^v I +' 1 ) i and O n ' test be the 
algebra generated by b # (f),f E h test , and l n . Note that O n ' test is 
dense in O n . Define 

Qtest ._ qZ ^ qKu test . . . QK n ,test ^ (gg) 

which is dense in O, and 

C:={att:aE O test } = T>(D~ r ), 

where D is the positive operator as defined in (80), section 4. We make 
the following additional assumption. 

(B4) The perturbation Hamiltonian V(s) E O test , for s € [0, 1]. 

In order to characterize the quasi-static evolution of nonequilibrium 
steady states, we introduce the new notion of an instantaneous NESS. 
Define an instantaneous NESS vector to be 

n g (s) := DU(-9)P Q g (s, 9)U(9)DVL . (84) 

Note that fl g from Corollary 4.4, section 4, has the same form as (84). 

It is important to note that introducing the operator D is needed 
to remove the complex deformation. 

We have the following Theorem, which effectively says that if a 
system, which is initially in a NESS, is perturbed slowly over a time 
scale t ^> tr, where tr is some generic time scale (tr = maxsgp,!] tr( s ), 
and Tft( s ) is the relaxation time to a NESS, see proof of Corollary 4.4 
in the Appendix), then the real state of the system is infinitesimally 
close to the instantaneous NESS, and the difference of the two states 
is bounded from above by a term of order 0(r _1 ). 
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Theorem 5.1 (Adiabatic Theorem for NESS) 

Suppose assumptions (Bl), (B3) and (B4) hold. Then there exists 
gi > 0, independent of s G [0,1], such that, for a G O test ,s G [0,1], 
and < \g\ < g±, the following estimate holds 

sup | <n ff (0), D _1 a t ™(a)n> - (ft fl (s), D _1 afi)| = 0(r _1 ) , (85) 
se[o,i] 

as r — > oo. 

Proof. Note that assumption (B3) implies assumption (B2), and 
hence the results of Theorem 4.3 about the spectrum of L*{t,9), for 
9 G I~(k) and fixed t G R, hold. The proof is now reduced to showing 
that the assumptions of Theorem 2.2 are satisfied. Choose 9 G I~(k). 
It follows from assumption (B3) and Lemma A.l in the Appendix, 
that the deformed C-Liouvillean L*(s,9) with common dense domain 

V = V(jCq) C\V(N) generates the propagator Ug T \s, s' ,9), s' < s, 
which is given by 

d s Uf\s,s',9) = -iTL* g (s,9)U^(s,s',9) , for a' < s;Uf\s,s,9) = 1. 

(86) 

This implies that (Al) and (A2) are satisfied. Furthermore, (A3) 
follows from the second resolvent identity 

(L* g (s,9)-z)- 1 = {U9)-z)-\l + gV to \s,9){U9)-z)- l r l , (87) 

and the results of Theorem 4.3, section 4. We also know that zero is an 
isolated simple eigenvalue of L*{s, 9) such that dist(0, a{L*(s, #))\{0}) > 
d, where d > is a constant independent of s G [0, 1]. This implies that 
assumption (A4) holds. Again using the resolvent equation (87) and 
assumption (B3), P^{s,Q) defined in (82) is twice differentiable as a 

bounded operator for all s G [0, 1], which imply (A5). Let Ua T \s, s', 9) 
(with domain V) be the propagator of the deformed adiabatic evolu- 
tion given by 

d s U£\s, s', 9) = -trL* a (s, 9)U^(s, s', 9) for s' < s ; u£\s, 8,6) = 1, 

(88) 

and 

Ll(s, 9) = i;(s, 9) + -[P g (s, 9), P g (s, 9)} . (89) 
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(Here, the () stands for differentiation with respect to s.) Since (Al)- 
(A5) are satisfied, the results of Theorem 2.2 hold. 

P°(s,6)U^( s ,0,6) = U^(s,0,9)P°(0,9) , (90) 

and 

snp \\U^(s,0,9)- [7^,0,0)11 = 0(0, (91) 
se[o,i] 

as r — > oo. 

For /i the single particle Hamiltonian of the free fermions, e lht 
leaves D(e k ^ p2+1 ) invariant. Therefore, for a G O test , a^a) G O tes ', 
where ckq corresponds to the free time evolution. Moreover, together 
with assumption (B4) and the boundedness of V, this implies (using 
a Dyson series expansion) that a T g s {a) G O test . 

Now, applying the time evolution on C(0,Q), and remembering 
that DQ = Q, U(9)Q = f2, the fact that U{9) and D commute, and 
the definition of the instantaneous NESS, it follows that 

(Q g (0), D'^ia)^ = (U^(s,0,0)P^0,9)n,a(9)Q) . (92) 

Using the results of Theorem 2.2, it follows that 

(Uf\s,0,0)P°(0,0na(9)n) 

= (U^(s, 0, 0)P°(O, 9)Q, a(9)Q) + 0(0 

= (P°(s,9)Ui T \s,0,9)n,a(9)n) +0(0 

= (p°(0^ r ^o,#)iMW + o(0 • 

The fact that (^ r) (s, 0, 9))*Q = Q implies 

DP°(s,9)Uf\s,0,9) = \n g (s,9))(n\Uf\s,0,9) 

= \n g ( s ,9))((uf\ s ,o,9)yn\ 
= \n g ( s ,9))(n\ = DP°(s,9). 

It follows that 

(n g (0),D- l a g TS \a)n) = (ng^D-iaty + Ofr- 1 ) , 
for large r. □ 
Remarks. 
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(1) Positivity of entropy production. If the interaction Hamiltonian 
gV(t) is time-periodic with period r, ie, V(t + r) = V(t), it 
is shown in [A-SF3] that the final state of the coupled system 
(introduced in section 3) converges to a time periodic state after 
very many periods. It is also shown that entropy production per 
cycle is strictly positive (Theorem 6.3 in [A-SF3]). The infinite 
period limit, r — > oo, is equivalent to the quasi-static limit. 
Hence, entropy production in the quasi-static evolution of NESS 
of the model considered in this paper is strictly positive. 

(2) An example of a reversible isothermal process. As a second appli- 
cation of Theorem 2.2 in quantum statistical mechanics, one may 
consider a concrete example of an isothermal process of a small 
system coupled to a single fermionic reservoir, and calculate an 
explicit rate of convergence (0(r -1 )) between the instantaneous 
equilibrium state and the true state of the system in the quasi- 
static limit t — > oo (see [A-SF1]). Under suitable assumptions 
on the form factors, one may show that there exists a constant 
<7i > such that, for a in a dense subset of O and < \g\ < g±, 
the following estimate holds 



as r — > oo, where p TS is the true state of the system at time 
t = rs, and oj@ s is the instantaneous equilibrium state, which 
corresponds to the zero eigenvalue of the time-dependent stan- 
dard Liouvillean. 

6 Appendix 

Existence of the deformed time evolution 

Choose 9 E I~(S), where 5 appears in assumption (B2), section 
3. The deformed time evolution is given by the propagator U 9 (t,t',9) 
which satisfies 




(93) 



d t U 9 (t,t',6) = -tL* g (t,6)U 9 (t,t',6) , U g (t,t,0) = 1 . 
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The following Lemma guarantees the existence of U g (t,t',9). Let 
V :=V(C )nV(N), 

and denote by 

C := sup sup \\V tot {t,9)\\ 
ten eel- (8) 



< 



V2 



sup £ |1 + e-&r 1/2 (Wi(t)\\H HS ,B) + \\e-^ /2 fi(t)\\H^,B)) < oc 



2 ten,zei(5) 



due to assumption (B2), section 3. 
Lemma A.l 

Assume (B2), choose 9 e I~ (5) U K and \g\ < g i} and fix i G R. 
Then 

(i) L*(t, 9) with domain!) generates a contraction semi-group e~ luL *^ t,6 \ a > 
on U. 

(ii) For ip G V, e~ iaL 3^il) is analytic in 9 E For 9' G R and 
9 g I- (6) UR, 

U(9')e~ i(rL ° m U(-9') = e - i ° L U t > 9 + e ') . 

(Hi) U g (t, t', 9)U g (t', t", 9) = U g (t, t " , 9) for t > t' > t" . 

(iv) U g (t,t' ,9)T> C V , and for ip G V, U g (t,t' ,9)ip is differentiable in 
t and t' such that 

d t U 9 (t, t', 9)4, = -iL* g (t, 9)U g (t, t', 9)4, , 

d t ,Ug(t,t',9)4, = iu g (t,t',e)L* g (t',9)i, . 

(v) For 9' G R ; 

U(9')U g (t, t', 9)U(-9') = U g (t, t', 9 + 9'). 
Moreover, U g (t,t',9) is analytic in 9 G I~(o~)- 
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Proof. Claim (i) follows from Phillip's Theorem for the pertur- 
bation of semigroups (see [Kal] chapter IX). Claim (ii) follows from 
assumption (B2), the resolvent identity 

(L* g (t, 9) -z)' l = (C (9) - z)-\\ + V to % 9)(C (9) - z)- 1 )- 1 , 

U(9')L* g (t,9)U(-9')=L* g (t,9 + 9'), 
and the fact that 

e -*L* g (t,9) = J_ f e -'*(i L *f t ,e) - Z)' l dz , 

2ni J r 

where T is a contour encircling the spectrum of L*(t, 9). 

Claims (iii) and (iv) are consequences of Kato's Theorem [Ka2], to 
which we refer the reader. Without loss of generality, rescale time such 
that t = ts,s E [0, 1], and let L* g n (sT,9) = L* g (r^9) for n E N\{0} 

and s E g, *±1], k = 0, • • • , n - 1. Moreover, define U£{ts, ts', 9) := 

e -ir(s-s')L^(r^e) j f J < g / < § < fc±l ? &nd {J^TS, TS 1 , 0) = E/£(rS, Tfi", #)1>(t 

if < s' < s" < s < 1. It follows from (ii) for 9' E R, that 
U(9')U^(rs, ts', 9)U{-9') = l>(rs, ts', 9 + 9') , 

and that U g (Ts,Ts' ,9) is analytic in 9 E where 5 appears in 

(B2). Claim (v) follows by taking the n — > oo limit (in norm). □ 



Glued Hilbert space representation 

We want to show that 

^(L 2 (R + ; B)) <g> F(L 2 (R + ; B)) = F(L 2 (R; B)) . 

Let Q be the vacuum state in the fermionic Fock space JF(L 2 (R + ; B)). 
For fermionic creation/annihilation operators on JF(L 2 (R + ; £>)), 

&#(/) : = y m{u)dudaf{u, a)b*(u, a) J E L 2 (R + ;B), 
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define the creation/annihilation operators on .F(L 2 (R + ; £>))©.7 : '(L 2 (R + ; £>)) 

as 

bf (/):=&#(/) ® 1 5 
&?(/) := (-1) N ®6 # (7), 

where 7 corresponds to complex conjugation. Note that 6/ and 6 r 
anti-commute. Let a and a* be the annihilation and creation oper- 
ators on the fermionic Fock space JF(L 2 (R+; £>) © L 2 (R + '; £>)) , such 
that they satisfy the usual CAR, and let Q be the vacuum state in 
J r (L 2 (R+; #) ©L 2 (R+; B)). An isomorphism between ^(L 2 (R+; B)) © 
F(L 2 (R + ; B)) and F(L 2 (R+; B) © L 2 (R+; B)) follows by the identifi- 
cation 

b*(f) = a*((f,0)), 
b*(g) = a#((0,g)), 

Now we claim that J r (L 2 (R + ;B) © L 2 (R + ;B)) is isomorphic to 
JF(L 2 (R; £>)). Consider the mapping 

j : L 2 (R+; B) © L 2 (R+; B) 9 (/,(?)—» /i € L 2 (R; B) , 

such that 

h(u, a) 



y/m(u)f(u,(r) ,u>0 
y/m(\u\)g(\u\,<7) ,u<0 



This mapping is an isometry, since 
Mf,9)\\hvL 



dudam(u)\f(u, a)\ 2 + / dudam(u)\g(u, a) 

J~R+;B 



R+,B Jn+;B 

= II/IIl 2 (r+ ; b) + IMIl 2 (r+;<b) • 

Moreover, the mapping j is an isomorphism, since, for given h G 
L 2 (R;B), there exists a mapping j" 1 : h — ► (f,g) e L 2 (R + ;B) © 
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L 2 (R + ;B), such that 



f(u,a) : 



1 



h(u, a), u > , 



g(u,a) : 




/i(H, a), u < . 



V m (\ u \) 



Proof of statements in Section 4 



Proof of Lemma 4-1 

C (9) restricted to the N — nl sector is 

4 n) (0) = C J] + s 1 + --- + s n + n9 , 



(94) 



where si, • • • , s„ are interpreted as one-particle multiplication opera- 
tors. For im# 7^ 0, it also follows from (94) that 



V = {4> = {^ (n) } : V W e T>(C { o\B)) and ^ ||4 n) (^)V (n) || 2 < oo}, 



and hence ^o(^) is a closed normal operator on P. Claims (ii) and (hi) 
follow from the corresponding statements on the sector N — nl. □ 

Proof of Lemma 4 ■ 2. 

The first claim (i) follows from the fact that gV tot (t, 9) is bounded 
for 9 G 1(5) due to assumption (B2) and the fact that the reservoirs are 
fermionic. It also follows from assumption (B2) that (g,0) — > L*(t,9) 
is analytic in 9 G Analyticity in g is obvious from (62). Assume 

that Re9 = 0. It follows from assumption (B2) that the resolvent 
formula 



(L^M)-^)- 1 = (Co(t,9)-z)- 1 (l+gV tot (t,9)(C (9)-zy 1 r l , (95) 



9 Although the results in this subsection are a very simple extension of those 
proven in [JP1,2,3] to the time-dependent case, they are sketched here so that 
the presentation is self-contained. The reader can refer to those references for 
additional details. 



n 
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holds for small g, as long as z belongs to the half-plane {z e C : < 
c < Imz}. Since (Co(t, 9) — z)^ 1 is uniformly bounded as ImO j for 
j£R and /mz large enough, and V tot (t,9) is bounded and analytic 
in 9, claim (iii) follows from the Neumann series expansion of the 
resolvent of L*(t,9). □ 

Proof of Theorem 4-3. (i) The resolvent formula 

(L;(t,9)-z)- 1 = {C {9)-z)-\l + gV to \t,9){U9)-z)- 1 )- 1 , (96) 

holds for small g and z in the half-plane {z e C : < c < Imz}. Note 
that 

II^MXA^)-*)- 1 !! < MCKA^)-^ 1 !! 



dist(z,ri(£ (9))) 



where C is given by (77) and i](£q(9)) is the closure of the numerical 
range of Cq. Fix #i such that it satisfies (76), and choose e such that 



e > ^fl > 0. Let 



fc-M 

G(z/, e) :={zGC: Imz; > v\ dist(z,rj(£ (9)) > e}. 



Then 



sup ||^(t,0)(£ o (0)-^~l<— , 



zeG(v,e) 9l 

uniformly in t. If \g\ < g±, the resolvent formula (96) holds on G(u, e), 
and, for m > 1, 



m— 1 



sup ||(z-L;(*,fl))- 1 -X;(«-A)(*,e))- 1 (»V tot (*,fl)( 2! -£o(e))- 1 ) , '|| < 

2£G(i/,€) i = Q 1 

(97) 

uniformly in t. It follows that 

|J G(v,e)cp(L* g (t,9)), (98) 
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where p(L*(t, 9)) is the resolvent set of L*(t, 9). Moreover, setting m = 
1 in (97), it follows that, for \g\ < gi/2, 

\\P g (t,9)-P \\<l, 

and hence P g (t, 9) is analytic in g. 

Fix (go,9o) G C x I~(5) such that \g \ < gi- Since L* o (t,9 ) 
and L* Q (t, 9) are unitarily equivalent if (9 — 9 ) e R and the discrete 
eigenvalues of L* o (t,9) are analytic functions with at most algebraic 
singularities in the neighbourhood of 9 , it follows that the pure point 
spectrum of L* Q (t, 9) is independent of 9. 

(ii) Analyticity of T g (t) in g follows directly from (i) and the definition 
of T g (t). Since ||T S (£) — 1|| < 1 for \g\ < gi/2, T g 1 ^ 2 (t) is also analytic 
in g. Inserting the Neumann series for the resolvent of L* {t, 9), gives 

oo 

T g {t) = l + Y.9 j T {j \t), (99) 
j'=i 

with 

T^\t) = j ^(«-£ E )- 1 Po^(t,fl)((«-£ (fi))- 1 V w («,fi)) , '- 1 Po(«-/: E )- 1 • 

(100) 

Similarly, 

oo 

M g (t) = £z + J2g J M^(t), (101) 

3=1 

with 

M ° )(t) = / ^- z ( z -^r lp oV tot (t, 9)((z-C (9))- 1 V tot (t, 9)) j - 1 P (z-C i: )- 1 . 

(102) 

The odd terms in the above two expansions are zero due to the fact 
that Pq projects onto the N = sector. The first non-trivial coefficient 
in the Taylor series of S 9 (t) is 

£( 2 )(t) = M^ 2 \t) - -(T {2 \t)C s + C*T {2 \t)) (103) 
= \ j t){Z ~ ^ + {Z ~ C ^~^ *)) ' < 104 ) 
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with 

£(z,t) = P V g tot (t,9)(z - c (e))- 1 v g to %9)p . 

□ 

Details of the calculation of the discrete spectrum ofL*(t,0) 

Denote by P k , k — 0, • • • ,3, the spectral projection onto the eigen- 
states of £ s , and let 

T k 2 \t):=P k ^ 2 \t)P k ,k = r -- ,3. 

Consider first the nondegenerate eigenvalues (E k = ^2uo,k = 2,3). 
Using the fact that 

1 1 
limine = VV-; 

e\0 X — it X 

lim/m = inSix), 

e\o x — ie 

and applying the Cauchy integration formula gives 

iter?) = -Yvvf d J JiMr \ 

3 2^ 7 R u-2uj q 

ImTf = -|£||£(2u*,f)|£, 

i 

and 

iter?) = --Tvvf d J UuMl , 
imr?) = —Y^WM^Wb ■ 

i 

Now apply degenerate perturbation theory for the zero eigenvalue. 
Using the definition of fi : p { and given in section 3, 

iter?] = o, 

TmV m _ n V \mu ,t)\\l ( e*-» -e*"»\ 
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Therefore, to second order in the coupling g, 

E 2t3 (g,t)= T (2u; - 1 -g 2 VV [ du—^—- J^UMWl) 

- i\? E ll^o,*)||| + 0(g 4 ), 

i 

while 

E 0>1 (g,t) =g 2 a ,i(t) + O(g 4 ) , 
where a ,i(t) are the eigenvalues of the matrix 

™ 2f 2cosh(/3 i w ) \-e~^ 

Since Q is an eigenvector corresponding to the isolated zero eigen- 
value of L g (t,9) (by construction, L g (t,9)Q = 0), then zero is also an 
isolated eigenvalue of L*{t,9). (One way of seeing this is to take the 
adjoint of the spectral projection of L g (t, 9) corresponding to Q, which 

is defined using contour integration.) Note that ip — ( \ ] is the eigen- 



vector corresponding to the zero eigenvalue of £( 2 )(t). Hence, 
E o (g,t)=0, 

Ek{g,t) = -tng 2 J2\\W^,t)\\l + 0(g 4 ) . 



Proof of Corollary 44 (NESS) 

Proof. Define k := min(^-, • • • , S), where S appears in assump- 
tion (B2), and let 9 G I~(k). We already know the spectrum of L*(t, 9) 
from Theorem 4.3. For a E O test , 

hm(n,a£(a)n) = lim (ft, e itL9 ae~ itLg Q) 

t— >oo y t—*oo 



= lim (e- itL *n, aVl) 

= lim (e- itL ^n, a(9)n) 

t— >oo 

= lim — ( / du(u + irj- vmy^-^+^n, a (e)n) , 
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for i] > 0. One may decompose the last integral into two parts (see for 
example [JP1]). The first part is 

lim —(ldz{z-L*{0))- l e~ izt n,a(9)n) = (Q g , D _1 aQ) , 
t^oo 2iri J 7 

where the zero-energy resonance is 

tt g := DU{-9)P° g {9)U{6)Dn = DU{-6)P° g {9)n . 

The second term converges to zero exponentially fast as t — > oo, since 

1 f°° 

— ( / {u-i{fi-€)-L* a (e))- 1 e- i(u - i ^- £))t n,a(e)n) = 0{e~^~ e,)t ) , 
where < e' < e < \Im0\ =: \i\ (see also Theorem 19.2 in [Rud]).D 
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